Magnetic effects on the phase transitions in 
unconventional superconductors 

Dimo I. Uzunov 

Max-Plank-Institut fiir Physik komplexer Systeme, Notnitzer Str. 38, 01187 Dresden, 

Germany, 

[Permanent address: CP Laboratory, Institute of Solid State Physics, Bulgarian 
Academy of Sciences, BG-1784 Sofia, Bulgaria.] 

Abstract 

The effect of magnetic fluctuations on the critical behaviour of unconventional 
ferromagnetic superconductors (UGe2, URhGe, etc.) and superfluids is investigated 
by the renormalization-group method. For the case of isotropic ferromagnetic order 
a new unusual critical behaviour is predicted. It is also shown that the uniaxial and 
bi-axial magnetic symmetries produce fluctuation driven phase transitions of first 
order. The results can be used in interpretations of experimental data and for a 
further development of the theory of critical phenomena in complex systems. 

Experiments ^1 |2] at low temperatures (T ~ 1 K) and high pressure (P ~ 1 GPa) 
demonstrated the existence of spin triplet superconducting states in the metallic com- 
pound UGe2. The superconductivity is triggered by the spontaneous magnetization of 
the ferromagnetic phase that occurs at much higher temperatures |Sl HI E]- The fer- 
romagnetic order coexists with the superconducting phase in the whole domain of its 
existence below T ~ 1.2. The same phenomenon of existence of superconductivity at low 
temperatures and high pressure in the domain of the (T, P) phase diagram where the 
ferromagnetic order is present was observed in URhGe jE]). These remarkable phenom- 
ena occur through phase transitions of first and second order and multi-critical points 
which present a considerable experimental and theoretical interest. A fragment of (-P, T) 
phase diagrams of itinerant ferromagnetic compounds jT] is sketched in fig. 1, where the 
lines Tp{P) and Tc{P) of the paramagnetic(P)-to-ferromagnetic(F) and ferromagnetic-to- 
coexistence phase(C) transitions are very close to each other and intersect at very low 
temperature or terminate at the absolute zero {Pq,0). At low temperature, where the 
phase transition lines are close enough to each other, the interaction between the real 
magnetization vector M{r) = {Mj{r);j = l,...,m} and the complex order parameter 
vector of the spin-triplet Cooper pairing [lOj, ipi^r) = {ipa{r) = {ip'^ + a = 1, ....n/2} 
{n = 6) cannot be neglected and, as shown here, this interaction produces new fluctu- 
ation phenomena. 

In this letter a new critical behavior for this type of systems is established and described. 
The new critical behaviour occurs in real systems with isotropic magnetic order but does 
not belong to any known universality class [7j. Thus it could be of considerable exper- 
imental and theoretical interest. Due to crystal and magnetic anisotropy a new type of 
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fluctuation-driven first order phase transitions occur, as shown in the present investiga- 
tion. The quantum effects [71 E] on the properties of these novel phase transitions are 
briefly discussed. 

Both thermal fluctuations at finite temperatures (T > 0) and quantum fluctuations (cor- 
relations) near the P-driven quantum phase transition at T = should be considered but 
at a first stage the quantum effects jH] can be neglected as irrelevant to finite tempera- 
ture phase transitions {Tp ~ Tc > 0). The present treatment of a recently derived free 
energy functional ^] by the standard Wilson-Fisher renormalization group (RG) [7j 
shows that unconventional ferromagnetic superconductors with an isotropic magnetic or- 
der (m = 3) exhibit a very special multi-critical behavior for any T > 0, whereas the 
magnetic anisotropy (m = 1,2) generates fluctuation-driven first order transitions [7]. 
Thus the phase transition properties of spin-triplet ferromagnetic superconductors are 
completely different from those predicted by mean field theories |3 EJ [TTl [T2I. The 
results can be used in the interpretation of experimental data for phase transitions in 
itinerant ferromagnetic compounds [T3] . 

The study presents for the first time an example of complex quantum criticality character- 
ized by a double-rate quantum critical dynamics. In the quantum limit (T 0) the fields 
M and ip have different dynamical exponents, zm and z^, and this leads to two different 
upper critical dimensions: dff = 6 — zm and = 6 — z^. The complete consideration of 
the quantum fluctuations of both fields M and ip requires a new RG approach in which 
one should either consider the difference {zm ~ z^) as an auxiliary small parameter or 
create a completely new theoretical paradigm of description. The considered problem is 
quite general and presents a challenge to the theory of quantum phase transitions [H]. The 
obtained results can be applied to any natural system with the same class of symmetry 
although this letter is based on the example of itinerant ferromagnetic compounds. 

The relevant part of the fluctuation Hamiltonian of unconventional ferromagnetic super- 
conductors in [m [12] can be written in the form 



where V ^ L'^ is the volume of the d— dimensional system, the length unit is chosen so that 
the wave vector k is confined below unity (0<A; = |fc| < l),5f>0isa coupling constant, 
describing the effect the scalar product of M and the vector product (ipxip*) for symmetry 
indices m = {n/2) = 3, and the parameters t {T — Tf) and r ~ (T — Tg) are expressed 
by the critical temperatures of the generic {g = 0) ferromagnetic and superconducting 
transitions. As mean field studies indicate [TT| lU]. Ts{P) is much lower than Tc{T) and 



The fourth order terms {M'^,\ip\'^, M'^\iIj\'^) in the total free energy (effective Hamilto- 
nian) [H [3 [HI [121 have not been included in eq. (1) as they are irrelevant to the present 
investigation. The simple dimensional analysis shows that the (yf— term in eq. (1) corre- 
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spends to a scaling factor b^-'^/'^ and, hence, becomes relevant below the upper borderline 
dimension = 6, while fourth order terms are scaled by a factor b'^~'^ as in the usual 
(/)^— theory and are relevant below < 4 (6 > 1 is a scaling number) [7j. Therefore we 
should perform the RG investigation in spatial dimensions d = 6 — e where the g-teim in 
eq. (1) describes the only relevant fluctuation interaction. Moreover, the total fluctuation 
Hamiltonian ^ E] contains off-diagonal terms of the form kikjipa4'p'i ^ 7^ J and/or 
a j3. Using a convenient loop expansion these terms can be completely integrated out 
from the partition function to show that they modify the parameters (r, t, g) of the theory 
but they do not affect the structure of the model (1). Such terms change auxiliary quan- 
tities, for example, the coordinates of the RG fixed points (FPs) but they do not affect 
the main RG results for the stability of the FPs and the values of the critical exponents. 
Here we ignore these off-diagonal terms. 

One may consider several cases: (i) uniaxial magnetic symmetry, M = (0,0,M3), (ii) 
tetragonal crystal symmetry when ip = [ipi, 1^2,0), (iii) magnetic order (Mi,M2,0), 
and (iv) the general case of cubic crystal symmetry and isotropic magnetic order (m = 
3) when all components of the three dimensional vectors M and ip may have nonzero 
equilibrium and fluctuation components. The latter case is of major interest to real 
systems where fluctuations of all components of the fields are possible despite the presence 
of spatial crystal and magnetic anisotropy that nullifies some of the equilibrium field 
components. In one-loop approximation, the RG analysis reveals different pictures for 
anisotropic (i)-(iii) and isotropic (iv) systems. As usual, a Gaussian ("trivial") FP {g* = 0) 
exists for all c? > and, as usual j7j, this FP is stable for d > 6 where the fluctuations 
are irrelevant. In the reminder of this letter the attention will be focussed on spatial 
dimensions d < 6, where the critical behavior is usually governed by nontrivial FPs 
{g* 7^ 0). In the cases (i)-(iii) only negative ( "unphysical" jTH]) FP values of g"^ have been 
obtained for d < 6. For example, in the case (i) the RG relation for g takes the form 

g' = b'-^/'-^>g{l+g^KMb), (2) 

where g' is the renormalized value of g, rj = {Kd-i/8)g'^ is the anomalous dimension 
(Fisher's exponent) [7j of the field M3; Ka = 2^''^n-'^/'^ /T{d/2). Using eq. (2) one obtains 
the FP coordinate [g"^)* = — QGvr^e. For d < 6 this FP is unphysical and does not describe 
any critical behavior. For d > 6 the same FP is physical but unstable towards the param- 
eter g as one may see from the positive value t/g = — lle/2 > of the respective stability 
exponent i/g defined by 6g' = b^^Sg. Therefore, a change of the order of the phase tran- 
sition from second order in mean-field approximation to a fluctuation-driven first order 
transition when the fluctuation (yf-interaction is taken into account, takes place. This con- 
clusion is supported by general concepts of RG theory |7] and by the particular property 
of these systems to exhibit first order phase transitions [H] in mean field approximation 
for broad variations of T and P. 

In the case (iv) of isotropic systems the RG equation for g is degenerate and the e- 
expansion breaks down. A similar situation is known from the theory of disordered 
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systems ^H] but here the physical mechanism and details of description are different. 
Namely for this degeneration one should consider the RG equations up to the two-loop 
order. The derivation of the two-loop terms in the RG equations is quite nontrivial because 
of the special symmetry properties of the interaction g-term in eq. (1). For example, some 
diagrams with opposite arrows of internal lines, as the couple shown in fig. (2), have oppo- 
site signs and compensate each other. The terms bringing contributions to the g-veitex 
are shown diagrammatically in fig. 3. The RG analysis is carried out by a completely 
new e^/'^-expansion for the FP values and e^^^-expansion for the critical exponents; again 
e = {6 — d). RG equations are quite lengthy and here only the equation for g is discussed. 
It has the form 

g' = 6(^-2r?^-^M)/2^ ^ j^g2 + 3(25 + C)g^] , (3) 

where 

A=^ [2\nb + e(ln6)2 + (1 - b^){2r + t)] , (4) 



B = {%\? - 1) - lllnfc - 6 (ln&)^] , C = [\nb + 2 {\nbf] , (5) 

rjM and i]^ are the anomalous dimensions of the fields M and ip, respectively. The one- 
loop approximation gives correct results to order e^^^ and the two-loop approximation 
brings such results up to order e. In eq. (4), r and t are small expansion quantities with 
equal FP values t* = r* = K^g'^. Using the condition for invariance of the two A;^-terms 
in eq. (1) one obtains rjM = Vi^ = where 

^ _ ^d-i o A 13 o 



Eq. (3) yields a new FP 



1 - ^^^-1^' ) • (6) 



g* = 8{37r'Y^\2e/13f\ (7) 

which corresponds to the critical exponent r] = 2(26/13)^/^ — 2e/3 (for d = 3, i] ^ —0.64). 

The eigenvalue problem for the RG stability matrix = [{dfii/d^j); {^i, ^2, fJ's) = (^^5^,5')] 
can be solved by the expansion of the matrix elements up to order e^^^. When the eigen- 
values \j = Aj{b)by^ of M. are calculated dangerous large terms of type 6^ and 6^(ln6), 
(6^1) 14j in the off-diagonal elements of the matrix A4 ensure the compensation of re- 
dundant large terms of the same type in the diagonal elements Aiu. This compensation is 
crucial for the validity of scaling for this type of critical behavior. Such a problem does not 
appear in standard cases of RG analysis [3 As in the usual 0^-theory ^1] the ampli- 
tudes A^- depend on the scaling factor 6: Ai = A2 = l+{27/13)b'^e, A3 = l-(81/52)e(ln6)2. 
The critical exponents yt = Vi, Ur = 2/2 and Ug = are 6-invariant: 

/2e 197 
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Vt = Vr — 18(2e/13)^/^, and = — e > for (i < 6. The correlation length critical 
exponents = l/Ur and vm = ^/Ut corresponding to the fields ip and M are 

1 5 /27 103 1 /27 5e . . 

These exponents describe a quite particular multi-critical behavior which differs from the 
numerous examples known so far. For c? = 3, z/^ = 0.78 which is somewhat above the usual 
value z/ ~ 0.6 -r- 0.7 near a standard phase transition of second order but um = 1-76 
at the same dimension {d = 3) is unusually large. The fact that the Fisher's exponent [7j 
1] is negative for d = 3 does not create troubles because such cases are known in complex 
systems, for example, in conventional superconductors The present e-expansion is 
valid under the conditions < 1, e^/^(ln6) <^ 1 provided 6 > 1. These conditions are 
stronger than those corresponding to the usual 0'^-theory [71 . This means that the 
present expansion in non-integer powers of e has a more restricted domain of validity than 
the standard e-expansion. Using the known relation [7j 7 = (2 — 77)1/, the susceptibility 
exponents for (i = 3 take the values 7^ = 2.06 and 7m = 4.65. These values exceed even 
those corresponding to the Hartree approximation |7^ (7 = 2z/ = 2 for d = 3) and can be 
easily distinguished in experiments. 

The critical behavior discussed so far may occur in a close vicinity of finite temperature 
multi-critical points (Tc = Tf > 0) in systems possessing the symmetry of the model (1). 
In certain systems, as shown in Fig. 1, this multi-critical points may occur at T = 0. 
In the quantum limit (T —>■ 0), or, more generally, in the low-temperature limit [T ^ 
fi] fi = [t,r);kB = 1] the thermal wavelengths of the fields M and exceed the inter- 
particle interaction radius and the quantum correlations fiuctuations become essential 
for the critical behavior [HI El- The quantum effects can be considered by RG analysis 
of a comprehensively generalized version of the model (1), namely, the action S of the 
referent quantum system. The generalized action is constructed with the help of the 
substitution (— 7i/T) S[M{q),ip{q)]. Now the description is given in terms of the 
(Bose) quantum fields M{q) and ip{q) which depend on the {d + l)-dimensional vector 
q = {ui, k); uji = 2-nlT is the Matsubara frequency (/i = 1; Z = 0, ±1, . . . ). The fc-sums 
in eq. (1) should be substituted by respective g-sums and the inverse bare correlation 
functions (r + k"^) and (t + k"^) in eq. (1) contain additional u;;— dependent terms, for 
example [SI IH| 

{\i^M)?r' = \^i\ + k^ + r. (10) 

The bare correlation function (|Mj(g)|)^ contains a term of type \uji\/k^ ^ where 6=1 and 
6 = 2 for clean and dirty itinerant ferromagnets, respectively [H]. The quantum dynamics 
of the field ip is described by the bare value 2; = 2 of the dynamical critical exponent 
z = whereas the quantum dynamics of the magnetization corresponds to zm = 3 (for 
= 1), or, to Zm = 4 (for 6 = 2). This means that the classical-to-quantum dimensional 
crossover at T — > is given hj d ^ {d + 2) and, hence, the system exhibits a simple mean 
field behavior for d > 4. Just below the upper quantum critical dimension d^^ = 4 the 
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relevant quantum effects at T = are represented by the field ip whereas the quantum 
[uji-) fluctuations of the magnetization are relevant for d < 3 (clean systems), or, for even 
for d < 2 (dirty limit) [Sj. This picture is confirmed by the analysis of singularities of the 
relevant perturbation integrals. Therefore, the quantum fluctuations of the field ip have 
a dominating role below spatial dimensions d < 4. 

Taking into account the quantum fluctuations of the field and completely neglecting the 
wj-dependence of the magnetization M, eo = (4 — (i)-analysis of the generalized action S 
has been performed within the one- loop approximation (order ej). In the classical limit 
(r/T <C 1) one re-derives the results already reported above together with an essentially 
new result, namely, the value of the dynamical exponent = 2 — (2e/13)^''^ which 
describes the quantum dynamics of the field ip. In the quantum limit [r/T ^ 1, T — > 
0) the static phase transition properties are affected by the quantum fluctuations, in 
particular, in isotropic systems {n/2 = m = 3). For this case, the one-loop RG equations 
corresponding to T = are not degenerate and give definite results. The RG equation 
for g, 



yields two FPs: (a) a Gaussian FP {g* = 0), which is unstable for d < 4, and (6) a FP 
(5'^)* = — 127r^eo which is unphysical [{g"^)* < 0] for (i < 4 and unstable for d> 4. Thus the 
new stable critical behavior corresponding to T > and d < 6 disappears in the quantum 
limit T — s> 0. At the absolute zero and any dimension d > the P— driven phase transition 
(Fig. 1) is of first order. This can be explained as a mere result of the limit T —* 0. The 
only role of the quantum effects is the creation of the new unphysical FP (fe). In fact, the 
referent classical system described by Ti. from eq. (1) also looses its stable FP (7) in the 
zero-temperature {classical) limit T ^ but does not generate any new FP because of 
the lack of g^-term in the equation for g'; see eq. (11). At T = the classical system has 
a purely mean field behavior P| which is characterized by a Gaussian FP {g* = 0) and is 
unstable towards T-perturbations for < (i < 6. This is a usual classical zero temperature 
behavior where the quantum correlations are ignored. For the standard 0^— theory this 
picture holds for d < 4. One may suppose that the quantum fluctuations of the field ip 
are not enough to ensure a stable quantum multi-critical behavior at Tc = Tp = and 
that the lack of such behavior is in result of neglecting the quantum fluctuations of M. 
One may try to take into account these quantum fluctuations by the special approaches 
from the theory of disordered systems, where additional expansion parameters are used 
to ensure the marginality of the fluctuating modes at the same borderline dimension du 
(see, e.g., Ref. [Hj). It may be conjectured that the techniques known from the theory of 
disordered systems with extended impurities cannot be straightforwardly apphed to the 
present problem and, perhaps, a completely new supposition should be introduced. 

In conclusion, the present results may be of use in interpretations of recent experi- 
ments PHI in UGe2, where the magnetic order is uniaxial (Ising symmetry) and the ex- 
perimental data, in accord with the present consideration, indicate that the C-P phase 




(11) 
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transition is of first order. Systems with isotropic magnetic order are needed for an 
experimental test of the new multi-critical behavior. 
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Figure 1: {P,T) diagram with a zero-temperature multicritical point {Po,0). Para- (P), 
ferromagnetic (F), and coexistence (C) phases, separated by the hnes Tf[P) and Tc{P) 
of P-F and F-C phase transitions, respectively. 




Figure 2: A sum of (^^-diagrams equal to zero. The thick and thin hnes correspond to 
correlation functions (iV'aP) and (|Mjp), respectively; vertices (•) represent gf-term in 
eq. (1). 




Figure 3: Diagrams for g' of third and fifth order in g. The arrows of the thick hnes have 
been omitted. 



